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of  the  Riccat.i  equation,  as  a  function  of  the  initial  condition.  , 

i  i 

liv  the  time  dependent  casoj,  the  domain  of  global  existence,  is  shown  to  Contain  a 
certain  cone.  > 


DD  :.™A473 


•  -a/m.'  w 

UNCLASSIFIED 
~~  "***  V»  »  t  la* 

i 


THE  GEOMETRY  OF  THE  RICCATI 
EQUATION 


by 

JOSE  MARIA  RODRIGUEZ-CANABAL 


A  Dissertation  Presented  to  the 
FACULTY  OF  THE  GRADUATE  SCHOOL 
UNIVERSITY  OF  SOUTHERN  CALIFORNIA 
In  Partial  Fulfillment  of  the 
Requirements  for  the  Degree 
DOCTOR  OF  PHILOSOPHY 


(Aerospace  Engineering) 


June  1972 

hF05/?  7 /-#rf/ 


#  Approved  f  r 

/&  Sistrifcuti  \\  u :l4l.  :• 


ACKNOWLEDGEMENTS 


The  author  wishes  to  egress  his  deep  gratitude  to  Professor 
R.  S.  Buoy  for  his  able  guid;  re,  inspiration  and  encouragement. 

Dr.  Bucy  has  been  very  gener  -ts  with  his  time  for  discussions  and 
teaching.  The  author  also  v  i..  ,ss  to  extend  his  sincere  thanks  to 
Professors  H.  K.  Cheng  and  A.  Troesch  for  their  interest. 

Thanks  are  also  due  to  Mr.  P.  D.  Weidman  who  provided  edito¬ 
rial  -issistance  and  to  the  author's  wife,  Eugenia,  for  her  typing. 

This  work  was  made  possible  by  the  support  of  the  Air  Force 
Office  of  Scientific  Research  under  grant  AFOSR-71-2141  and  the 
support  of  the  Army  Research  Office,  Durham  under  grant  DA-ARO-D- 
31-12**172"G*»8. 


i  t 


Approved  for  put  lie  release- 
distribution  unlimited. 


CONTENTS 


Chapter 

I  INTRODUCTION 

II  FIXED  POINTS 

III  GEOMETRY  OF  THE  RICCATI  EQUATION 

IV  2-D  RICCATI  EQUATION 


Page 

I 

5 

19 

29 


REFERENCES 


33 


NOTATION  CONVENTION 


A 

|Aj 

MI 

I 

Re  A 

{«!•••) 

a-*>b 

cm 

n 

€ 

4 

c 


dA 

dt 

absolute  value  of  A 
appropriate  norm  of  A 
identity  matrix 
real  part  of  A  (p.  15) 
set  of  a's 

a  tends  to  b,  or  a  changes  into  b 

(2  x  2)  matrix 

set  intersection 

empty  set 

belongs  to 

does  not  belong  to 

set  inclusion 


iv 


CHAPTER  I 


INTRODUCTION 

The  Importance  of  the  Rlccatl  equation  (RE)  in  itself  and 
in  connection  with  the  optimal  linear-quadratic  control  and  Kalman- 
Bucy  filtering  theory,  has  motivated  several  works  examining  the 
behavior  of  its  solution.  In  particular,  results  pertaining  to 
sufficient  conditions  for  global  existence,  monotoneity  properties, 
asymptotic  character,  invariant  sets  and  equilibrium  points  have 
been  reported  in  the  literature.  A  discussion  of  the  interesting 
problem  of  a  priori  bounds  for  the  solution  is  presented  in  (4J. 

As  Willems  l 1 9}  has  pointed  out,  various  aspects  of  the  RE 
need  further  investigation.  Specifically,  one  is  interested  in 
i)  performance  criterion  which  are  not  positive,  ii)  the  cost  linear 
functional  in  the  control,  iii)  control  with  conflicting  aims,  and 
iv)  conjugate  point  structure  in  variational  problems. 

This  latter  case  is  of  particular  importance  in  the 
numerical  solution  of  the  RE:  although  the  solution  of  the  differ¬ 
ence  matrix  equation  may  be  well  behaved,  in  the  sense  of  being 
finite  for  all  time  and  convergent  to  an  equilibrium  solution,  it 
does  not  pertain  to  the  problem  under  consideration  since  the 
solution  lies  outside  the  interval  of  di sconjugacy.  The  only  indica¬ 
tion  that  one  gets  from  the  numerical  computation  is  a  change  of 


1 


2 


the  definiteness  of  the  solution  between  two  consecutives  iterations. 

There  exists  a  duality  principle  between  the  optimal  linear- 
quadratic  control  problem  and  the  optimal  l inear-gaussian  filtering 
problem.  The  duality  is  achieved  by  a  simple  change  of  the  dynamics, 
sensor,  and  actuator  along  with  a  reversal  of  time.  See  [10, 

Section  5 ]  .  This  correspondence  allows  one  to  formulate  results 
in  either  control  or  filtering  form  3'd  obtain  results  for  the 
dual  problem  by  a  reinterpretation  within  the  appropriate  context. 

Henceforth  we  will  be  concerned  with  solutions  of  the  RE 
with  time  running  forward  and  the  initial  condition  imposed  at 
the  starting  time.  We  will  also  concentrate  on  the  autonomous  case, 
i.e.,  when  the  coefficient  matrices  are  time  independent.  For  this 
case,  under  certain  regularity  conditions,  it  has  been  shown  [6] 
that  there  exists  a  fixed  point  of  the  RE,  called  P_,  which 
generates  a  convex  cone  in  the  partial  ordering  induced  by  the 
positive  semi-definite  matrices,  such  that  if  the  initial  condition 
belongs  to  this  cone  then  the  solution  exists  globally.  The  result 
was  proved  by  the  application  of  a  theorem  of  Reid  [16J. 

The  importance  of  P_  in  connection  with  the  structure  of 
the  solution  of  the  RE  was  also  recognized  in  (19)  whereas  its 
physical  interpretation  as  the  steady  state  error  variance  of  the 
present  state  given  future  observations  was  made  clear  in  [4). 

Also  a  solution  [4,  Corollary  3 • 2 ]  of  the  RE  is  given  when  the 
initial  condition  tends  to  infinity;  this  is  a  particular  case  of 


i 


the  general  form  presented  in  Chapter  ill. 

t 

The  main  results  determined  in  this  analysis  are:  > 

i)  Necessary  and  sufficient  conditions  for  global  existence 
of  the  RE  as  a  function  of  the  initial  condition  for  the 

!  5 

autonomous  case. 

i 

ii)  Sufficient  conditions  for  the  solution  to  be  in. the 

domain  of  global  existence  for  the  time  dependent  case. 

iii)  A  general  existence  theorem  for  the  equilibrium  points 

! 

P  and  p  .  (The  reader  is  referred  to  Chapter  II  for  the 
definition  of  P+l  P_  and  Pq) . 

iv)  A  method  for  finding  the  equilibria. 

There  are  several  techniques  available  for  the  computations 

l 

of  the  points  P+f  P_  and  Pq .  One  vay  to  find  P+  is  by  the  method 
of  quasi-linearization  or  Newton's  method  in  function  space.  What 

i 

results  is  a  numerical  computation  [5,  Chapter  8],  [12]  which  can 
efficiently  compute  P+  ;  P_  can  also  be  determined  by  iterating 
to  a  related  equation  l  6].  This  method  is  not  applicable  for  the 
determination  of  Py .  , 

A  second  method  rakes  use  of  information  contained  in  the 
Hamiltonian  matrix,}/,  associated  with  the  RE.  One  form  of  this 
technique,  due  to  Potter  [l*»],  is  valid  when}/  is  diagonal  izable;. 
Martensson  [l3i  has  extended  the  method  to  include  a  general  , 

Hamiltonian  matrix.  In  the  present  work  we  extend  another  form  of 


i 
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4 

the  Hamiltonian  method,  due  to  Bass  and  Roth,  [l],  {  53  >  and  |58j, 
which  is  based  on  a  partition  of  'the  characteristic  polynomial  of  K 

i 

The  use  of  the  Hamiltonian  provides  not  only  P,  and  P  ,  but 

T  " 

also  any  , ,  if  it  exists,  as  well,  as  results  concerning  the  number 
of  equilibria.  'In  [ 1 9 j  a  geometric  method  is  presented  which, 
however,  is  not  amenable  to  computation. 

’  The  wor;k  on  determination  and  structure  of  equilibria 

1 .  > 

presented  here  was  done  independently  of  that  by  Willems  [l9j  (see 

comment  in  [6j).!  We  note  that  although  he  has  presented  some  results 

; 

i 

concerning  the  equilibria  of  the  RE,  many  of  these  assertions  are 
without  proof.  Fu/thermore,  $e  feel  our  proofs  are  simple  and 

i  t 

direc't.  In  addition  our  method  is  numerically  tractable  because  it 
is  eigenvector  free,  whereas’ those  of  [  1 3  J ,  114)",  and  [19]  are 
numerically  impractical.  In  fact  we  have  a  digital  program  which 
effectively  computes  all  real  symmetric  equilibria  and  is  quite 
accurate. 


*  Note  that  (14)  concerns  the  determination  of  complex  equilibria. 


CHAPTER  II 


FIXED  POINTS 

For  real  symmetric  matrices  we  will  use  positive  definite, 
A>0,  positive  semi-def i ni te,  A^O,  and  strictly  positive  semi- 
definite,  A^O,  to  denote  that  the  matrix  A  has  all  its  eigenvalues 
greater  than  zero,  greater  or  equal  to  zero,  and  greater  or  equal  to 
zero  but  at  least  one  equal  to  zero,  respectively.  The  obvious  dual 
definitions  and  notations  are  used  for  negative  definite,  semi- 
definite,  and  strictly  semi-def ini te.  Nonnegative  definite  will  mean 
a  real  symmetric  matrix  which  is  either  positive  semi-definite  or 
indefinite,  rather  than  the  normal  use  with  the  meaning  of  positive 
semi -def i ni te  in  the  present  form. 

The  set  of  positive  semi -def i ni te  matrices  is  a  convex  cone,  ^ 
which  induces  the  following  partial  ordering  over  the  space  of  real 
symmetric  matrices: 

A  >  B  iff  A-B£0,  i  .e.  A  -  Bfi  & 

A  maximal  (minimal)  element  of  the  set  B  is  an  element  Pe  B  , 
such  that  there  is  no  element  in  B  -fp)  which  is  greater  (smaller) 
than  P.  Finally,  the  supremum  (infimum)  element  of  the  set  B  is  a 
P«B  such  that  P  is  greater  (smaller)  than  any  other  clement  of  B>  . 

Notice  that  in  the  space  of  symmetric  matrices  it  is  not 

possible  to  define  the  supremum  as  the  least  upper  bound,  since  the 

5 
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common  upper  bounds  of  two  given  matrices  may  be  not  comparable 
among  them,  for  instance,  given 


A  = 

0  2 
2  0 

B  - 

o  o' 
0  0. 

c  & 

‘i  r 

I  l 

D  = 

i  i 
Ll  2j 

then  C^A,  C^B  and  D^A,  D^B  but  C  and  D  are  not  comparable. See  [ 2 j . 


If  A  is  a  finite  real  symmetric  matrix  we  define  the  pseudo- 
supremum  of  A  end  0,  the  zero  matrix,  as: 

ps  sup  (A , 0}  =  T'  I.  A;  v  0]  T 

where  T'AT  =  diag  [^;]  ,(V'0Ja  diag[max  U;,  0) ]  ,  and  then 
ps  sup(A.B)  =  ps  sup(B,A)  =  ps  sup(A-B.O)  +  B 

0 

As  is  normal  3 &  and  &  will  denote  the  boundary  and  the 
interior  of  the  set  &  ,  respectively. 


We  consider  now  the  autonomous  RE: 

3Y  =  PF‘  +  FP  -  ph’hp  +  a 
P(t0)  -  r 


(2.1) 


where  prime  denotes  transpose  and  wich  Mj  n(R)  the  space  of  (1  x  n) 
real  matrices,  FeMnjn(R),  H6MS)n(R),  Q  =Q'eMn>n(R)  ,  T=r'e  Mnn(R) , 

P  =  P'e  Hn,n(R)-  The  right  hand  side  of  (2.1)  will  be  denoted  by 
S  (P) ,  and  let  TT(t,r)  be  the  locally  unique  solution  of  (2.1),  see 
[5,  Chapter  5j- 


The  Hamiltonian  matrix  associated  with  (2.1)  is 


h'h 

F 


(2.2) 


For  the  definition  of  uniform  asympiotic  stability  (u.a.s.) 


see  {5>Chapter  ! 3 • 
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A  simple  substitution  gives  the  following: 

I  DENT ITY  2.1.  For  X+X - 1  ,  X  a  real  number  and  P  j ,  P2  any  two  real 
symmetric  matrices 

S  *  A  Pj  +  A  P2)  =*S  (P,)  +  A  S(P2)  +X%P,,  P2) 
where  R(Pj,  PQ  >.  0  for  every  Pj,  P2  • 

DEFINITION  2.1. 

s+  «*  {r=r|  S ( r )  >  0} 
s_  ={r=r'i  s  cr)  ^  0} 

Notice  that  they  are  extensions  of  S+  and  S_  defined  in  (5>  Chapter  5) 
Here  we  do  not  make  the  restriction  that  Fe  &  . 

REMARK  2.1.  S+  is  convex.  It  follows  from  Identity  2.1. 

In  [19]  the  convexity  is  proved  by  the  use  of  the  linear  matrix 
inequality  (LMl)  and  a  complicated  string  of  equivalences.  A  direct 
proof  from  the  LMl  can  be  constructed  by  using  results  iri  f 9 ) . 

REMARK  2.2.  If  P,e  S+ ,  P S_  (P2e  Sj  andX<0  then  P=XPj+XP2e  S_  , 

(P  e  S_)  where  \  =  1 . 

THEOREM  2.1.  j_f  Ffi  S+l  (Fe S_)  then  JT( t ,  F  ,  to)  a  solution  of  (2.1) 
belongs  to  S+  ,  (S_)  for  al 1  time  tel  (interval  of  existence)  and 
is  monotone  increasing  (decreasing), 

PROOF:  See  L 5 .  Chapter  5]« 

THEOREM  2.2.  J_f  Tf(  t ,  f7 ,  t0 )  a  solution  of  (2.1)  exists  in  the  interval 


3  then  for  £  >  17  >  IT  ( t ,  F  ,  tc )  exists  in  the  interval  3  and 


TT ( t , ,t„)  >T[(t,f7  ,  t„)  for  te  'i 

PROOF :  See  1 16,  Lemma  2.3]. 
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DEFINITION  2.2.  p 
iff  S  (P)  =  0. 


p  6  Mn>r,(R)  is  an  ecJu*  librium  point  of  (2.1) 


R.EHARK  2.3.^  If  {Pi  fs  the  set  of  equilibrium  points  of  (2.1)  then 
{P)£S+n  S_. 

— °REM  2,3 Pl  *  (p2)  is,  a  maximal  minimal  element  of  S+  then 
P1  ,^P2^  -‘s  an  equilibrium  point  of  (2.1). 

PROOF!  Let  P,  be  a  maximal  element  of  S+  ,  then  by  Theorem  2.1., 
valid  with  Q.  =  Q.  ,TT(t  ,P))  >  P]  ,  hence  TT(t,Pj)  =  P)  forft:>  0  by 
maximal ity  and  therefore  P, «  {  P)  .  For  P£  the  proof  follows  by 
reversing  the  time  i.e.  t-*--  t. 

The  following  three  lemmas  ere  basic  for  the  proof  of  the 

main  theorem  concerning  the  existence  of  maximal  and  minimal  elements 

of  S  . 

+ 

— EMMA  —  ^  &  be_.a  closed  convex  cone  with  vertex  at  zero,  and 

S4-  bejajclos^  !_  convex  set  such  that  OeS^.  J_f  S  0  &  does  not  have  a 

maximal  element  in  the  partial  ordering  induced  bv  &  then  there 

exits  an  Ae  &  such.  j.  ha'  AAfiS+n&  for  every  real  number  a  >.  n 

PR00Fj_  Since  S +n  £  does  not  have  a  maximal  element  it  is  possible 

to  find  a  totally  ordered  sequence  An  with  no  finite  upper  bound, 

such  that  Ane  S+n  & ,  for  otherwise  by  Zorn's  Lemma  S  n£  has  a 

+ 

maximal  element.  Without  loss  of  generality  An  may  be  taken 
increasing. 
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Since  0 e § .  ,  there  exists  a  sphere  of  radius  k  about  0, 

'•  i 

j 

&k  =  {A  |  ||AIU  k}  £  S+  and^Bk/2n&  is  compact.  Define  Xn  so  that 

Cn  =  >An  An  and  II  Cn  II  =  k/2  .  '  ’ 

By  compactness  ofdBj^^  &  there  exists  a  subsequence  Cm  which 

converges  to  C  i.e.  Cm-»  c<!Vh 

Let  Dm  =  Cm+Am  <€  -  Cm)  with  /\»l  and  II  Dn'  -  C  1 1  -6  S  k/2 

©  ’>  } 

with  £  fixed.  By  construction  Dme  S+  as  Dm  is  in  the  €-sphere  about  ■ 

C  which  lies  in  the  k-sphere  about  0. 

i 

Consider  the  two  one-parameter  families  1 

L(Pm)  -  0„  +  '  0J  for  OS^I  .  ’ 

1 

M(  )  -AC  for  A  ^  0  . 

I 

L(Pm)es+  for  Vm,  since  S+  is!convex  and  Dm  Ame  S+  .  ‘ 

Suppose  there  exists  a  /\j  ,  such  that  M(  ^j)f  ,  and 

|  J 

let  choose  such  that: 


then  L(PJ  =/<m(1  -PJ  C  6  S+0  ^  ,  and  since  SpO^-l)  I  i  C  -  Cm  1 1 
then  00  .  » 

*  1  I  ! 

Suppose  there  exists  a  subsequence  /tm,  such  that 
/imA]  ■/6m'^kl<co  ’  this  lmPUcs  that  ^mr'  >  and  fi n:),Amrco 
contradicting  the  choice  of  , 

4 

Hence  there  exists  m0  large  enough, 'fqr  which/^  (1  -  /®m)  > ^  i  - 
and  ^C  e  S+o  £?  by  convexity  of  S+  ,  contradicting  the  supposition 
that  M(Aj)^S+n£?  ,  therefore  "XCe  S+()&  for  V/UO. 


LEMMA  2.2.  A  bounded  monotone  sequence  converges. 


PROOF:  See  t 173- 


I 


I 


1  (  *  •  '  ■  ?0 
LEMMA  2.3-  I f  S+  associated^  wi th  (2.1),  has;  a  maximal  element  P j 

ft  o 

!  then  there  exists  a  matrix  ?2>  Pj  such  that  S_  .  : 

PROOF:  Since  Pj  is  a  maximal  element  of  S+  then  S(P)i^  for  vF^p  . 
Let  1a  =  Pj  +  A  with  A^O  then:  ‘  <  1 

S(fi  j  =  S*  (A)  =  AFj  +  F,A  AH'HA^  £  for  ,V  A>0,  where 

I  i 

F)  =  F  -  PiH'H  .  • 

*  *  *1 
\  { 
t  Consider  the  perturbed  operator  ,  > 

s? (A) »  ~~  =  AFj  +  F.A  :  AH'HA'  +  kl  (2.3) 

M  1  1 0  ^  •  1 

with  k  a  positive  real  number  and  I  the'  (n  x  n)  identity-matrix. 

-  1  ' 

I  o  I  *  0  , 

S+  corresponding  to  (2.3)  is  not  empty  as  A  =  0  £  S+  . 

1  :  :  1 

Suppose  that'S.Ag’  associated  with  (2.3)  does  not  have  a 

max'imal  element,  thep  by  Lemma  2.1.  there  exists  an  A j  >  0  such  that 
,S*(*  A,')->  0  forV^O,  but  ,  ’  , 

,  S*UA)  «  >-S*  (A)  +  X(1  -A  )AH'HA  +  kl 
wi  th  S*  (A)^  £?'.  Let  ^mjn  =  niin[X)  A  eigegvblue  of  sf  (A))  <  0 
hei^ce  for  "Xj  >  k/l  Xmjnl  we  have  that  S*(X]A)£S+  ,  contradiction, 
therefore  the  set  S+0 (3  associated  with  (2.3)  has  a  maximal  element  ' 

•  •  1  i 

A  j  >  0  and  since  (2.3)  is  completely  control  1 abl e 'we  can  apply 
[5.  Theorem  i>.l)  so  that  if  JT( t , 0, t0 ) ■ i s  the  solution  of' (2.3) 

1  '  r-i*  ''  * 

with  ini t ( a  I  cond i t ion  0,  then  Tl ( t , 0, t0 ) >  0  for  t  >  t0  and  by 

I  ^ 

Theorem  2.2.  Aj  >TT  ( c  ,0,  t0)  therefore.  Aj  >  0,  and  since  TT(t,0,to) 

1  *  h 

is  a  bounded  monotone  sequence  it  converges  to  A  such  th'at  Aj£  A  >0 

and  sl(A‘)  =  ,0,  hence  P^'  -  Pj  +  A* 

l  * 

,  The  next  theorem  shows  that  there  is  at  most  one  maximal 
(minimal)  element  Pj  ,  of  S+  and  that 


!  1 


s+c{r=r'|r<P)}  n  {r=r'jr>P  j 


( minimal ) 


THEOREM  2.4.  I f  P|  ,  (Pj)  is  a  max  i  maP'el  erne  nt  of  S^.  then  it  is 
the  supremum  (infimum)  element  of  S+  . 

PROOF :  Suppose  BeS+  and  B  is  not  comparable  v.*i  th  Pj  .  By  Lemma  2. 3. 
there  exists  A,>P)  and  A}e  .  Then  with  “X +\  «■  l,*<0,  we  have 


XPj  +  X  Aj>  Pj  and  for 


X 


^max(B  "Aj)/^mjn(Aj"Pj)  ,^j<0, 


then  A]Pj  +  “Xj  A  |  >  B  and  by  Remark  2.2. .  P^  =  +  ^jAj  e  S  . 

■  TT(t,  Pj.to)  is  monotone  decreasing  as  t  increases,  Theorem  2.1, 
and  bounded  below  by  Pj  ,  Theorem  2.2.,  hence  converges  to  an  equi¬ 
librium  point  which  has  to  be  P)  ,  since  otherwise  Pj  will  not  be 
a  maximal  element  of  S+  . 

n(t.B  , t0)  is  monotone  increasing  and  bounded  above  by 
J7(t,P3,t  )  for  al 1  time,  hence  Pj>P£  >  contradicting  the 
assumption  that  they  were  incomparable.  The  proof  for  the  minimal 
element  goes  simmilarly  hy  reversing  the  time. 


Once  the  supremum  and  infimum  of  S^.  are  r.haracter i zed  as 
equilibrium  points,  when  they  exist,  will  be  denoted  by  P+  ,  P_  , 
respectively.  The  use  of  this  clasical  notation  will  be  justified 
latter. 


We  now  give  necessary  and  sufficient  conditions  for  the 
existence  of  P+  and  of  the  P  . 

THEOREM  2.5.  P+  ,  ( P_ )  exists  iff  S+  is  non-empty  and  there  exists 
a  matrix  Lfs;  (Lfs-)such  that  F  -  LjH'H  ,  (-F  +  L2H'H)  is  u.a.s. 
PROOF^  Suppose  LeS.  and  F^  =  F  -  LH'h  is  u.a.s.  then  for  V  K e $+ 
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FjA  +AF[  =  AH*HA  +  A  (2.4) 

with  h=  K  -  L  and  A>0. 

Suppose  A0  a  solution  of  (2.4)  and  A0is  not  negative 
semi-definite,  then 

FlA+AF^  =  AH‘HA+  A  (2.5) 

has  by  assumption  at  least  one  solution  A0,  but  the  solution  to 
(2.5)  is  unique  as  F  is  u.a.s.  and  negative  semi-definite  since 

“*>_  -t 

A  =  -  eFLt  [aoH'HA0+  A  ]  eFlt  dt 

•'o 

Contradiction.  Therefore  A<0  and  L>:KforVKeS,  i.e.  L  is  an 
upper  bound  for  S+  which  implies  the  existence  of  a  maximal  element 
of  S  by  Zorn's  Lemma  and  by  Theorem  2.4.  of  P  . 

T  *T 

If  there  exists  P  ,  then  by  Lemma  2.3  there  exists  and 

T 

L>  P,  and  such  that  Le  S  ,  and  it  is  easy  to  see  that 

T  " 

Fl  (L  -  P+)  +  (L  -  P+)  F^  =  S(L)  -  (P+-L)  H'H(P+  -  L)  <  0 
and  Fl  is  u.a.s.  by  Liapounov's  Theorem. 

For  P_  the  proof  goes  simmilarly  and  is  omitted. 

The  following  two  lemmas  show  why  the  notation  P+  and  P_ 
has  been  used  although  no  definiteness  for  them  has  been  assumed. 
LEMMA  2.4.  F  -  P.h'h,  -F  +  P  h'h  are  u.a.s,  iff  S.£  0 

r  '  "  *T*  “  *T 

PROOF :  It  follows  from  Liapounov's  Theorem.  See[l9»  Theorem  5j • 
LEMMA  2.5-  With  S.^-  <i>  ,  P  ,  and  P  are  the  only  equilibrium 

'  I  --  r  I  T  '  *T*  mmm^mmm  mm  '  '  .  _ 

solution  of  (2.1)  that  make  F ,  =  F  -P.H'H  ,  -F  =  -F+  P  H'H  u.a.s. 

T  T  “  * 

PROOF:  Suppose  K^-P.  makes  F„  u.a.s.,  then 

+  I\ 
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F+(K-P+)  +  (K-P+)F+  -  (K-P+)H'H(K-P+)  =  0 
implies  K  -  P+<0.  (See  proof  of  Theorem  2.5.),  and 

Fk(P+-K)  +  (P+-K)F^  -  (P^tiH'H^-K)  =  0 
implies  P  -  K<  0,  therefore  K  =  P  .  See  also  [19,  Lemma  2j  for 

T  T 

another  proof. 

REMARK  2.U.  Theorem  2 A  or  the  proof  of  Lemma  2.5  shows  that  if 
P6  is  any  other  equilibrium  point  different  from  P+  id  P_,  then 
P  <  Pe  <  P+;  in  Chapter  III  we  wi  1 1  show  that  P_  ^  Pe  ^  P+. 

REMARK  2.5.  In  the  case  of  Q.  =  GG1 ,  i .e. ,  Q  >  0,  then  P  =  0  €  S 
and  hence  if  there  exists  P+  it  has  to  be  positive  semi -def i ni te. 
With  that  constrain  the  following  results  are  known. 

i)  If  (F,H)  is  completely  observable  and  (F,G)  completely 
controllable,  then  there  exists  an  unique  positive  semi- 
definite  equilibrium  point;  it  is  positive  definite  and 
is  P.,  and  the  matrix  F  -  P  H 1 H  is  u.a.s.  See  [  5],  I  1 0 J . 

ii)  If  (F.H)  is  completely  observable,  then  1 im  TT(t,0,to) 
exists  as  t-*  <»  and  equal  to  P.  Further  S(P)  =  0  and 
P  is  positive  semi-def i ni te.  See  [5]>  [10).  For  Q  =  Q1, 
[19]. 

ill)  If  (F h')  is  stabilUable  andfG1,  F1)  is  observable, 
then  lim  7T  ( t ,  0 ,  t0 )  =  P+  as  t-*  <»  and  P+  is  the  unique 
positive  semi-definite  solution  of  S(P)  =  0.  Further  P+ 
is  positive  definite  and  the  matrix  F  -  P+H'H  is  u.a.s. 
See  [20]. 


iv)  If  (F1,  H1)  is  stabilizable  and  (G  1 ,  F1)  Is  detectable, 
then  S(P)  =  0  has  at  least  one  solution  P  positive 
semi-def ini  te  and  the  matrix  F  -  PH'H  is  u.a.s.  See  f 20] . 
v)  If  (F1,  H1)  is  stabilizable,  then  TT(t,0, t0)  is  bounded 
on  [t0,oo). 


Notice  that  Theorem  2.5  includes  the  above  results  as 
particular  cases  but  is  more  general  since  no  assumption  on  the 
definiteness  of  Q.  has  been  made,  as  in  the  following  example. 


EXAMPLE  2.1.  Let 


F  = 


■1  0 
1  -l 


H  =  (1  ,  0)  Q  = 


-1  1 
1  -3 


S+  is  non-empty  as  P  =  fg  _2|  e  s+>  L  =  0  e  S_  and  FL  is  u.a.s. , 

therefore  there  exists  P .  ,  in  fact  P.  -  ,°r 

+ ’  +  0-1.5 

unbounded  below  since  there  is  no  l.  e  S+  such  that  -F[_  is  u.a.s 
Indeed  S4  =  {  ~]Q  °  |  c  «  -1.5  ) 


.  However  S+  is 


The  next  problem  to  be  considered  is  how  to  effectively 
compute  the  equilibrium  points  of  (2.1).  Notice  that  any  method 
depending  on  eigenvectors  is  numerically  troublesome  as  is  well 
known  from  perturbation  theory,  see  [ll,  Chapter  l]  . 


We  will  use  fol lowing  [ 1  ]  ,  l  ]  the  Hami 1 tonian  matrix  . 
LEMMA  2.6.  det(  1{  -Xl)  =  R(‘A)  =  (-l)nA(A)  A  (-A) 
where  A(A)  is  a  polynomial  of  degree  n. 

PROOF:  See  f  5, Chapter  8  )  . 


DEFINITION  2.3-  A  factorization  of  R(A)  is  any  set  (X| ,  An) 
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of  eigenvalues  of  °H ,  possibly  repeated,  such  that  R{^)  = 

(-1)  A  (X)  A  (*A)  where  A  (X)  is  a  polynomial  of  degree  n  with 
roots  (^\j  ,^2*  •  •  *^n)  • 

THEOREM  2.6.  Every  equilibrium  point  P  of  (2.1)  satisfies: 

(-P  ,  I)A(W)  =  0  (2.6) 


where  AP.)  corresponds  to  some  factorization  and  A  (^\)  = 
det(F  -  PH'H  -  >1) . 

PROOF:  If  P  is  an  equilibrium  point  of  (2.1),  let  T  =  [*? 


I 

P 


then 


f’X  T 


F  -  PH'H  0 

h’h  -f‘  +  H'HP 


Let  A  P)  be  the  characteristic  polynomial  of  F  -  Ph'h  which 
corresponds  to  some  factorization  of  R(^\),  then  A  (70  =  ‘ 
and 


A  B 

T  = 

-PB  +  D  C  -  PA  +  (D  -  PB) P 

C  D 

B  A  +  BP 

A(p  -  PH’H)  0 

*  A(-F  +  H'HP) 

therefore  D  -  PB  and  C  =  PA  and  the  theorem  follows.  See  [5, Theorem 

8. 5  J  . 

COROLLARY  2.1.  P+,  (  P_  )  is  given  by  (2.6)  when  the  factorization 

is  such  that  A  U)  is  a  Hurwitz,  (totally  non-Hcrwi  tz)  polynomial , 
i  .e. ,  Re  \  j  <  0,  (  Re  A;  >  0  )  for  V  i ,  where  A  p;)  =0 


It  is  now  possible  to  know  how  many  equilibria  may  exist. 
If  all  eigenvalues  are  distinct  and  real,  from  each  it  is 
possible  to  form  two  different  factorizations,  namely  one  with  X 
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the  other  with  -  "X  ,  then  there  are  2n  different  factorizations  and 
therefore  there  are  2n  equilibrium  points,  since  the  half  lower 
part  of  the  eigenvectors  of  expand  Rn. 


EXAMPLE  2.2. 

F  -  0  1  H’H  -  1  0  0  -  1  0 

p  [o  !J  H  H  “  [o  ij  ^  "  [0  3. 

then  RO)  =  -  5  A2  +  4  ,  and  (-P  ,  1)  Aj(K)  =  0,  i  =  1,2,3,** 

gives: 


1  0 

-1  0 

)  o' 

-1  0 

0  3 

P_  = 

0  -1 

Pe, 

0-1 

pe2  " 

0  3 

If  all  eigenvalues  are  real  and  one  of  them  is  repeated, 

say  of  order  s,  then  the  (n  -  s)  different  ones  can  be  combined 
n_$ 

in  2  ways  and  the  s  repeated  gives  (s  +  l)  possible  combinations, 
hence  either  there  are  (s  +  l)2n  s  equilibria  as  in  example  2.3. 

If  however  the  minimal  polynomial  of  coincides  with  one  A. (  ) 
corresponding  to  some  factorization,  then  (2.6)  is  satisfied  for 
every  P,  and  yet  P  has  to  satisfy  A  (F  -  PH'm)  =  0  and  there  is 
the  posibility  of  parametric  families  of  solutions  as  in  Example  2.4. 


EXAMPLE  2.3. 


h'  H  =  Q,  = 


I  I 
1  1 


From  the  different  factorizations  we  get 
P+  -  1  P-  "  -I  Pe 


1  0 
0  -l 


However  R(A)  -  (^2  -  l)2,  since  ^ 


l  is  not  a  minimal  polynomial 


of  7{  • 


EXAMPLE  2,4. 


H'H  =  Q  - 


l  0 
0  1 
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In  this  case  P+  =  I  and  P„  =  -]  .  For  Pe  ,  Ae(X)  is  the  .minimal  i 
polynomial  of  K  ,  although  the  only  P  that  satisfy  Ag(F  -  PH'h)  =  0 

,  a  family  of  equilibria. 


are  Pe  = 


cose  sine 
sine  -cos e 


The  following  example  shows  that  P+.  P„  may  be  singular 
and  yet  P+  -  P_  be  positive  definite,  i.e.,  F+  is  u.a.s. 
EXAMPLE  2.5. 

“  h  nl  .  h  nl 

Q.  » 


F  « 


then  P+  = 


2  0 
0  0 


I  0 
0  -1 


H  H  = 


and  P 


[0  o' 

-  -  [o  -2 


1  0 
0  1 


and  P+  -  P_>  0 


If  there  is  a  complex  eigenvalue  and  the  others  real -and 
distinct,  then  each  factori2ation  has  to  include  that  eigenvalue 
and  its  complex  conjugate  as  A  (A)  is  the  characteristic  .polynomial 
of  a  real  matrix.  Then  there  are  2n  '  different  factorizations. 


EXAMPLE  2.6. 


Inthis  case  there  are  cniy  P+  s  I  and  f>_  =  -1 

If  finally,  there  is  a  complex  repeated  eigenvalue  of  crdef 
s  and  the  c.hers  real  and  distinct,  it  is  possible  to  have 
(s  +  l)2n  different  factorization  and  posibilities  isolated 
eui librium  points  or  families  of  solutions. 


The  above  consideration  are  simpl  ookeeping  of  the 
combinatorial  ways  to  factor  the  r  f  iyenvalues. 

In  the  general  case  of  real  '  s  5 1  ict,real  and  repeated, 
complex  distinct  and  complex  repeated,  we  can  count  the  possible 
number  of  i  motorizations  by  considering  the  different  classes. 


i 


flEMAR'K  2.6.  If,  [-P  ,  l]  A  (X)  =  0  has  no  solution  for  a  certain 
factorization,  then  there  is  no  equilibrium  point  P  of  S(P)  =  0 
such'that  A(F  -  PH*H)  =  0,  since  (2.6)  is  in  general  a  necessary 
condition. 


.  The  following  are  examples  of  those  cases. 

! 

EXAMPLE  2.7'.  ■ 

1  o]  .  f  1  0|  0  0 

p  =  h'h  =  a  = 

l-J  -IJ  [o  0j  [p  1. 

the  eigenvalues  are  +1  ,!-1,  repeated  and  for 

r  2  -i] 


A  0\)  ,=  +2^  +  1  one  gets  ,P.  = 


-1  ’I 


A  (^)  =  “X  -  l'  one  gets,  through  the  use  of  the  characteristic 


polynomial,  of  Fe 


:  [6 

pe  =• 


e  i 


e 

0  +e)2 


For  A(^)  =  'X  -  2  V  +  1  ,  (2.6)  is  a  set  of  incompatible  equations 
and  there  is  no  ,P_.  Indeed  Theorem  2.5  asserted  it  and  S+  is 
unbounded  below.  . 


EXAMPLE  2.8. 


H  H  = 


The  set  of  equations  (2.6)  is  i.ncompat ible  for  the  different  factoriza¬ 
tions  and  therefore  there  are  no  equilibrium  points.  It  is  interesting 
to  note  that  S+  is  non-empty  but  S_  =  <j>  . 
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CHAPTER  ill 


GEOMETRY  OF  THE  RICCAT!  EQUATION 

.In  the  following  we  will  consider  mainly  the  solution  of  the 
RE,  in  the  form  (2.2)  of  the  previous  chapter,  i.e., 

=  PF1  +  FP  -  PH'HP  +  Q. 

(3-D 

P(t0)  =  r  =  r  e  Mn>n(R) 

In  this  case  it  is  possible  to  find  necessary  and  sufficient 
conditions  on  P  for  the  global  existence  of  TT(t,P).  The  importance 
of  P  in  connection  with  this  problem,  it  has  gone  unnoticed  until 
very  recently  1 6 ] ,  f  19],  maybe  due  to  the  fact  that  the  principal 
tool  used  was  a  solution  connected  with  the  Hamiltonian  matrix 
associated  with  (3-0  and  due  to  Radon  \l5]- 

ASSUMPTION  Al.  S.  is  bounded  above  and  below 

DEFINITION  3.1. 

j?p_  =  {r=r'l  r-p.  >  0) 

ep_  =  jr=r'|  r-P.  >°) 

^p_=(r.r'|  r-R^o) 

simi larly, 

\  ,(4)=fr=r|  r‘Pt-  (c,°} 

Finally  let  (Fe  )0£(<  be  the  set  of  all  other  equilibrium  points  of 
(3.1)  different  from  P  and  P  and  this  set  can  be  empty,  see 

T  “ 
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example  2.6. 

DEFINITION  3.2. 

Se| r-  r'|Te„n_n(R)j 

is  an  invarian  set  under  the  transformation  defined  by  the  RE  if 
for  Fe  B  then  P(t,r)£  B  f°r  V  t  £0. 


ASSUMPTION  A2.  P+  -  P_  is  non-singular. 

LEMMA  3.1.  Suppose  A1  and  A2  hold  and  F -P_  is  non-singuiar  then 
the  solution  of  ( 3 ■  1 )  i  s  given  for  t  e  "J  (the  interval  of 
existence  ofH(t.r))  by 

P(0»P.-«{(P+-PJ’1+ f‘,(t,te)[(r-PJ’1-(P+-Pj",]yr"'(t.t  r1)  1 

(3.2) 

where  (t,t0)  is  the  fundamental  matrix  of  F_=F-P-H  H 


PROOF;  The  proposed  solution  satisfies  (3-1)  since  defferentat ion 


gives 

p=(p-p_){f!  f"[(r-pj“1-(p+-pj“1]t  ,+f"l( 

F_](P-P_) 


or 

P=(P-P_){F_'[(P-P_)"1-(P+-P_)"1J+[(P-P_)"1-(P+-P_)"1jF.j(P-Pj 
where  a  rearranged  form  of  (3.1)  has  been  used,  namely 
P=(P-P_)  F'_  +  F_  (P-Pj  -  (P-P_)H'H(P-P_) 

and  since  P  (t,P+)  =  0  for  Vt  then  (3.2)  satisfies  (3-1).  The 
initial  condition  P  (t0)  =Fis  satisfied  as  ^(t0,t0)  =  ] 


THEOREM  3- 1 •  Under  the  conditions  of  Lemma  3-1 
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i )  JjF  F €  8p  the  solution  FI  (  t  ,  P  )  exists  globally  and  for 

t-ro  , JT(t, r )  —  p+ 

i i )  if  r  f  Sp  and  P -P_  is  non-singular,  then  the  solution 

TT(t,r)  has  a  finite  escape  time. 

PROOF :  To  prove  i)  we  have  to  show  that  P(t,T  )  is  finite  for  all 
time  t  or  equivalently,  that 

B(t)  -  (p+-p_)">f'/L(r-p_)"l-(p..-p_)_,]f'' 

is  non-singular  for  te  (to,00). 

Lemma  2 .h  gives  the  u.a.s.  condition  of  -F_  hence  (t,t0)^0 
as  t  —  m 
Let 

c(0  =  (p+-pJ_1-  f\ P+“P_Pf  ' 

rhen  C(t)  is  the  solution  of 
C  =  -F'  C  •  C  F_  +  H'H 

c(t.)  =  0 

Therefore 

0(t)  =  J*  H'H  f  \t,s)  ds  >  0 

and 

B ( t )  -  -P_)"'jr*,+C(t)>  0  as  r-P_>  0 

as  t-»-°o  ^  '-*•  0  ,  and  P(t,r)-,-P+ 

To  prove  part  ii)  we  note  that  B(t0)  is  not  positive  definite 
or  semi-definite  and  for  t-«-co  ,  B(t)  limits  to  a  positive  definite 
matrix,  namely  (P  -P  )  hence  as  the  eigenvalues  are  continuous 

T  “ 

functions  of  t,  for  some  finite  value  t  =  T,  B(T)  will  be  singular 


and  the  solution  will  escape  to  infinity. 


We  have  investigated  the  behavior  of  the  solution  for 
r  =r,e  n(R)  and  P  ~P_  non-singular.  To  study  the  case  of 
singular  F  -P_  we  need  some  more  results. 

LEMMA  3.2.  (Dua  I  of  Lemma  3.1.).  Suppose  A1  and  A2  hold  and  p 
non-singular  then  the  solution  of  (3-1)  Is  given  for  tel  ,  by 

p(t)=p++|  $''(t,tc)[  (r-p+)'1+(p+-p_)_1J  t0)-(p+-p_)"1  J 


;re  (j5(t,t0)  is  the  fundamental  matrix  of  F-P^'H 


PROOF:  Same  as  Lemma  3-1. 


This  Lemma  generalizes  Corollary  3.2  of  [4]. 


LEMMA  3.3.  (Complementary)  Let  ([AH  =A-P 


A 1 ,  A2  hold  and 


V  S  o  £p  then  the  solution  of  (3-1)  is  given  for  tel  by 

p(t)=p_+  f  irjfyVjf ]f-  1[p+3  _1)  trFj-'Srj^f' 

(3.4) 

wi  th  'f'  and  J  as  before. 

For  r€{F  =  rir-P.goJ  the  solution  is  given  by(3.4)  replacing 

trl*  b*  (-p#. 

PROOF:  Essentially  same  method  as  in  Lemma  3-1  with  the  use  of  the 
positive  square  root  of  positive  semi-def i ni te  matrices. 

REMARK  3.1.  Similarly  there  is  a  dual  complementary  lemma  for 
T-P,  being  singular. 


i 
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>  i 

LEMMA  3-4.  Suppose  A 1 ,  A2  hold  then  ( Pe )  £ ^  &p_ 

O 

PROOF:  For  T e  Tr  ,  P(t,F)-*P+  as  t-t®  hence  there  cannot  be  > 

o  :  1 

equilibria  other  than  P  in  Cp_  .  • 

+  i 

For  F 4  Sp.  and  T-P  ^  non-singular,  P(t,T)  escapes  to  ihfinity 
in  finite  time,  therefore  there  cannot  be  equilibria  .in  this 
region.  We  are  left  then  with  the  set  of  F  s  such  that  I  -P_  is 
singular.  Consider  the  subset  T -P_^  0  of  that  >  callit  SR  6p_ 

i 

(boundary  of  the  rear  cone  with  vextex  at  P_) .  Since  by  the  , 
change  t -*■  -t  the  structure  of  the  equilibrium  points  does  not 

•  ‘  i  1 

change  f 8] ,  however,  the  u.a.s.  ones  change  to  totally  non-u. a. s. 
and  vice-versa,  then  by  considering  the  solution  of  this,new  < 
equation,  o  K  up.  is  inside  of  the  region  of  convergence  of  the 

i 

solution  to  the  point  P  .  Therefore  there  is  not  equilibria  on 
‘d  ?v  £p_  ,  hence  (P@)  C  c) 

I 

LEMMA  3  ■  5.  Suppose  A 1 ,  A2  hold ,  then  ( Pp  ) c  ^  } 

!  i 

PROOF:  Note  that  the  transformation  t-*-  -t  makes  this  lemma  the 
dual  of  the  previous  one. 

l 

1  1  '  ,  1 
THEOREM  3.2.  Suppose  A1 ,  A2  hold,  then  (Pa)c  ^  ^  fl  d  Cp, 

PROOF:  Follows  from  Lemma  3-^*  and  3.5-.  ,  i  , 

REMARK  3-2.  I  f  Pe  6  (Pe)  then  P_$  PQ  ^  P+ 

^  t  5 

There  are  several  invariant  sets,  e.g  S.  and  S  .  Also 

1  “  X 

since  A  being  a  singular  positive  semi-definite  matrix  implies  A2 

is  singular,  then  from  Lemma  3-3  and  its  dual,  it  fo.llows  that 


1 


I 


! 


{r  =T  |  (T  -P_)  singular] 


{r=r'  |  p-.p  )  singular] 


are  invariant  sets. 

:  »i  i 

THEOREM  3-3-  :The  set  "3  S’ p_  H  1)  is  invariant  under  the  , 

I  1  i 

■  transformation  defined  by  the  RE.  .  , 

PROOF Follows  from  the  invariance  of  the,  above, sets. 

? 

I  i 

,  We  can  sumarize  the  previous  result  in  a  graphic  way.  Refer' 

,  '  '  '  ,  .  ' 

to  Fig.  3. I ,  where  the  structure  of  a  2  x  2  RE  is  represented  via 

the  space  R  with  axis  P  ^ ,  P  P^.  Fig.  3-1  is  just  a  suggstion, 

l 

'  since  actually  the  cones  are  cones  with  a  maximum  angle  of  TT/2  at, 

'  i  i  ,  i 

the  vertex  and  elliptical  cross  section  with  excentricity  e=l/  V27 

i  1  !  . 

,  We  can  define  the  domain  of  attraction  of  a  setB,  hs  the 

■  i 

set  of  F's  such  that  P(t,r)e  B  as  t-~«‘  .  Then  the  domain  of 

‘  O  1  ' 

attraction  of  P+  is  c\>_  and  the  one  of  7)  H  3  ffP+  js  3  {?p_ 

We  see  that  if  P  ,  P  ,  (P'  -P  )  '  exist,  the  signature  of  0, 

T  **  T  "  i 

and  P  are  immaterial  for  the  behavior  of  the  solution. 

‘  '  ’  ‘  . 

REMARK  3.3.  The  previous  results  are  a  remarkable  generalization 
of  the  one-dimensional  case, 

2  2  i 

p  =  2fp-hp+q  ,  '  '  ' 

where  the  domain  of  attraction  of  p+  is  (p  ,co)  as  depicted  in 
Fig.  3.2. 


i 


FIG.  3.2 


We  can  now  better  understand  some  of  the  unusual  features 
which  appear  in  the  numerical  solution  of  the  RE.,  see[l3]. 

if  r*  net  ,T  )  will  converge  to  P+  as  t  gets  large,  in 

spite  of  round-off  errors. 

If  T  ^  £?p_  the  solution  will  escape  to  infinity  in  a 
finite  time. 

If  F  e  ^  ,  theoretical  lyTT(t,n)e<)6R  for  Vtand 

tends  to  fi  9  as  t  •■*•«»;  this  is  the  only  way  of  reaching  a 
Pe  and  thus  provides  a  cure  or  in  some  cases  even  an  analytic 
solution.  See  chapter  IV.  Excep  for  very  special  cases,  the  round¬ 
off  error  will  throw  the  solution  outside  "d  £p_  ,  and  will  therefore 
escape  in  finite  time  or  converge  to  P  . 

This  explains  why  in  [13]  the  numerical  solution  of  (2.1) 
is  unable  to  reach  Pe  with  initial  values  in  the  domain  of 
attraction  of  Pe  . 

We  now  consider  a  time  dependent  RE. 
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p(t)  =  P ( t ) F 1 (t)  +  F(t)P(t)  -  P(t)H'(t)H(t)P(t)  +  d(t) 

p(t0)  =r  =  r'€  Mn>n(R) 

wi  th  Q.(t)  =  Q.  (t) . 


(3.5) 


By  using  the  fundamental  matrix  of  F(t)  it  is  possible  to 
change  (3.5)  into: 

Pj (t)  =  -p,(t)  Hj(t)  h, (t)  pj(t)  +  a,(t) 

Pi(t.)  -  r 

where 

W(t, t0)  =jto$(s,t)H'  (s)H(s)  <j'(s,t)ds  =jto  Hj(s)Hj  (s)ds 
C(t,t0)  =J^$(t,s)  a(s)  <^(t,s)ds  a  (s)ds 

with  W(t,t0),  C(t,t0)  the  observability  and  controllability  matrix, 
respectively.  See  [3]  . 

Hence  without  loss  of  generality  we  will  consider  a  RE 


wi th  F(t)  =0 

P(t)  =  -P(t)  M(t)  P(t)  +  Q(t) 

P(t„)  =r  =r'e  Mn>n(R) 

and  M(t)>°»  Q  =  Q*  continuous  functions  of  time,  and  let: 

M"  =  ps  sup  M(t) 
telt  o,  co> 

Q*  =  ps  inf  Q.(t) 

te[ta,co) 


(3.6) 


Notice  that  ps  sup  and  ps  inf  can  be  replaced  by  any  other 
convinient  upper  and  lower  bounds,  respectively. 

Consider  the  following  constant  RE: 

p* .  -p*  h"  .  <e  0.7) 

p-(t.)  =  r 


and  assume  ther<*  exists  the  minimum  element  of  S+  associated  with 


it,  V,  . 
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THEOREM  3.^.  I  f  M(t)  >0,  Q(t)  =  Q.'  (t)  are  continuous  functions  of  t 
and  there  exists  ^  ,  then  fl(tfr,t0),  solution  of  (3.6),  exists 

J_n  te[t0,oo)  for  T  >  H 

» 

PROOF:  t0)  solution  of  (3.7)  exists  in  te[t0,co),  from 

Theorem  3.1  and  with  Tj  =  Tl(t ,  T  ,  t0 )  -  71  (t.f  ,t0)  ,  then  t|  satisfies 

^  =  -  rf MTj  -  t|mtt'"'  -  rjnij  +  (a  -  q*)  -  n"(M  -  m")  rr 

ij(t.)  =  0 

therefore  Yj  (t)  >0  for  all  t,  i  .e. ,  f[(t,  P  ,t6)  >  (t,  r  ,ta) . 


CHAPTER  IV 


2-0  R I C CAT  I  EQUATION 

Through  this  research  numerical  computations  and  theoretical 
results  for  the  2-dimens lonal  RE  were  useful  since  the  space  M2  2(R) 
Is  Isomorphic  with  R^  over  the  field  of  real  numbers  and  addition  of 
matrices  and  vectors,  for  symmetric  matrices. 

The  equation  of  the  motion  of  the  solution  of  the  RE  with 
Te  3  0^  wj  i  ]  presented  here.  It  is  remarkable  the  natural  way 

in  which  improper  rotations,  [7],  appear  in  the  two-dimensional 
case.  It  should  be  expected  that  something  similar  will  happen  in 
a  higner  dimension,  although  we  did  not  investigate  this  case  any 
further. 

I  DENT  I T Y  k,  1 ,  For  (2x2)  re?l  matrices  the  following  identities 
arc  va 1 i d; 

i)  det(A±B)  B  det(A)  +  dct(B)  t  (trace  A) (trace  B)+tracc(AB) 
or,  _[_f  A  is  non-singular: 

ii)  det(AtB)  B  dct(A)  +  det(B)  -  dct(A)  tracc(A~l  B) 
where  dci(C)  B  determinant  of  C. 

LEMMA  ^i.i.  J_f  P+  -  P_  Is  non-singular  then  for  (2  x  2)  real  sym- 
rnetrlc  matrices  the  solut  ion  of: 
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Is  given  by; 


P«=  P.  +■ 


det  (P  -  P+)  =  0 
det  (P  -  P+)  =  0 

(1  +  Te) 


trace  {(P+  -  P^f1"  (I  +  Te  )} 

cose  sine 
sine  -cose 

:  Let  Pe  “  P-  =|b  c]  *  sat!sfy  ^e  second  equation  of  (4.1) 


-(4.1) 


(4.2) 


where: 


PROOF 


T  = 


and  P0  -  P_>0  the  following  relations  have  to  be  met: 

2 

ac  ■  b  =0,  a,  c^O  (4.3) 

This  family  of  solutions  can  be  represented  by  the  two-parameter 
function: 

Pe  -  P„  "Ml  +  Te) 

where  ^  is  a  positive  real  number  and  0e[O,  21X]  . 

Indeed  the  transformation: 

a  =  ^  ( 1  +  cos  6  ) 
b  =  ^  sin  0 
c  =  X  ( 1  -  cos  9  ) 

satisfies  relations  (4,3)  and  (a,  b,  c)-*-(^\,0)  is  one-to-one 

a  +  c 
*=  — 

2  b  \ 


0  =  s  in 


-1 


a  +  c  i 

where  0  is  fully  determined  from  the  other  relations,  and  is  well 
2b' 

'  £1  . 


defined  since 


a  +  c 


Now  by  the  use  of  Identity  4.1.  ii),  it  is  possible  to 
determine  the  value  of  X  so  that  the  first  equation  of  (4.1)  is 
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also  satisfied  and  the  proof  is  complete. 

,  I 

I 

MOTION  0N2Sp.fl9gp.: 

Sustitution  of  (?».2)  into  the  RE  gives:  ; 

0  0  Te  ©(I  +  Te)  trace {(P+  -  P.)“^jTeJ 

trace{(P+  -  Pj^U  +  Te)j  [trace  {(P+  -  K)’]  U  +  Te)}]2 

(k.k) 

F-  (I  +  Te)  +  (I  +  Te)  F-  (1  +  Te  )h'h(1  +  Te )  '  > 

trace{(P+  -  Pj'!(I  +  Te )}  [trace{(P+  -  P-.)”*1  (I  +  Te)}]2 

i 

where 

j_ r °  >i 

°  [-1  o 

Since  7)^036^  is  an  invariant  set  under  the  transformation 
defined  by  the  RE,  the  three  component  equations  derived  from  (^.*0 
should  be  non-independent;  therefore  we  are  permitted  to  choose  the 

i 

most  convinient  linear  combination  for  the  solution. 


CASE  I:  trace  }(P+  -  Pj'ljV}-*  0  *  , 

By  taking  the  trace  of  (*4.4)  the  solution  of  the  RE  with 
initial  condition  £p+  is  given  by:  ■ 


P(t)  =  P„ 


(1  +  Te ) 

trace  {(P+  -  Pj"1  (I  +  T* )} 


i 


where  8  sat i sf ies: 

^  _-trace{F_0  +  Te )}  trace{(P+  -  P_)_,(I  +  Te)}+  trace{(l  +  Te )H 'h] 

trace{(P+  -  P.)_IJ' Te] 

0(to)  =  So 


where  trace(AB)  =  trace(BA)  and  trace(j'Te)  =  0  ,  have  been  used. 


i 


i 


CASE  II:  trace{(P+' -  P.)“,j'Te}=  0  * 

i  • 

In  this  case  (4. It)  reduces  to  , 

j  ‘  ’  i  ■  (I  +  T0 )  H  F 

■91  «  F_(I  +  Te  )Te  J  +  (1  +  Ta)F.  TeJ 


trace  {(P+  -  F 

or  taking  the  trace  and  since  TeJ'Te  =  J  ,  then  the 

I  ; 

as  in  case  I  with  9  satisfying: 

I  ’  ’  ’  ♦ 

9  =!  trace;{F_(I  +  Te)  J} 
i  0  (tfl )  =  So 


i 
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(1  +  Te)Tej' 

_rT(I  +  Te)} 

solution  is  given 


I 
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